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Abstract 

A multifractal-like representation for multi-time multi-scale veloc- 
ity correlation in turbulence and dynamical turbulent models is pro- 
posed. The importance of subleading contributions to time correla- 
tions is highlighted. The fulfillment of the dynamical constraints due 
to the equations of motion is thoroughly discussed. The prediction 
stemming from this representation are tested within the framework of 
shell models for turbulence. 
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1 Introduction 



Turbulent flows are characterized by a highly chaotic and intermittent trans- 
fer of fluctuations from the stirring length, outer scale, Lq, down to the 
viscous dissipation length, inner scale. Id- The Reynolds number defines the 
ratio between the outer and the inner scales: Lo/ld = Re^^^. We talk about 
fully developed turbulent flows in the limit Re — > oo, in this limit it is safely 
assumed that there exists an inertial range of scales, la ^ r Lq, where the 
time evolution feels only the non-linear terms of the Navier-Stokes eqs. 

The highly chaotic and intermittent transfer of energy leads to non-trivial 
correlation among fluctuations of the velocity fields at different scales and at 
different time-delays 0] . 

The natural set of observable which one would like to control are the 
following: 

C^'^ir,R\t) = {6vm-5vU0)) (1) 

where 6vrit) = v{x + r, t) —v{x, t) and -C r < i? -C Lq. In (|l|) we have, for 
sake of simplicity, neglected the vectorial and tensorial dependencies in the 
velocity fields and velocity correlations respectively. 

Some subclasses of the Multi-Scale Multi-Time (MSMT) correlation func- 
tions (0) have recently attracted the attention of many scientists [0, |, |, |]. 
By evaluating (P with r = R, at changing R, and at zero-time delay, t = 0, 
we have the celebrated Structure Functions (SF) of order q + p. Further, we 
may also investigate Multi-Scale (MS) correlation functions when we have 
different lengths involved r ^ R at zero delay, t = as well as single-scale 
correlation functions (CF) by fixing r = i? at varying time delay t etc... 

Structure functions have been, so far, the most studied turbulent quanti- 
ties (see for a recent theoretical and experimental review). On the other 
hand, only recently some theoretical and experimental efforts have been done 
in order to understand the time properties of single scale (CF) correlations, 
(jp,q (j-^j- g] the scaling properties of multi-scale correlations (MS) 
at zero time delay, 0^'" (r, |0) 0, |, §. 

In this Paper, we propose and check a general phenomenological frame- 
work capable to capture all the above mentioned correlation functions and 
in agreement with the typical structure of non-linear terms of Navier-Stokes 
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eqs. In Section ^ the framework of the multifractal description of correla- 
tions is briefly recalled and critically examined. In Section ^ a representation 
for single scale time correlations is introduced, and its predictions are tested 
within the framework of shell models of turbulence. In Section ^ we deal 
with the most generic two-scales time correlation. 

2 Background: the multifractal description 
of time correlations 

One of the most important outcomes of experimental and theoretical anal- 
ysis of turbulent flows is the spectacular ability of simple multifractal phe- 
nomenology |]I|, to capture the leading behavior of structure functions and 
of multi-scale correlation functions at zero-time delays p|, ^ 0. This may 
appear not surprising because, as far as time-delays are not concerned, one 
may expect that (many) different phenomenological descriptions may well 
reproduce scaling laws typical of (SF) and of (MS) functions: multifractals 
being just one of these descriptions. More striking were the recent findings 
that multifractal phenomenology may easily be extended to the time-domain 
such as to give a precise prediction on the behavior of the time properties of 
single scale correlations. As soon as time enters in the game, one must ask 
consistency with the equation of motion: the major break-through was that 
one may write a time-multifractal description in agreement with the equation 
of motion^. 

Let us now quickly enter in the details of previous findings in order to 
clarify both the phenomenological framework and the notation that we will 
use in the following. 

Let us remind that the multifractal (Parisi-Frisch) description of single- 
time correlation functions is based on the assumptions that inertial range 
statistics is fully determined by a cascade process conditioned to some large 
scale configuration: 

6vr = Wir, R) ■ 6vR (2) 

"'^An important remark is now in order: when we refer to time-properties of turbulent 
flows we always mean the time-properties of the velocity fields once the trivial sweeping 
effects of large scale on small scales is removed, for example by choosing to work in a 
quasi-Lagrangian reference frame |^ . 
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where the fluctuating function W{r, R) can be expressed in terms of a super- 
position of local scaling solution W{r,R) ~ {^)^^^^ with a scaling exponent 
h{x) which assumes different values h in a class of interwoven fractal sets 
with fractal codimension Z{h) = 3 — D{h). From this assumption one can 
write the expression for any structure functions of order m, which in our 
notation [m = p + q) becomes: 

S'^iR) = C^^'^ {R, R\0) ~ {W {R, LoD {U^) (3) 

where we have introduced the shorthand notation d^jij^^ih) = dh Ij^j to 
define the probability of having a local exponent h connecting fluctuations 
between scales R and Lq. In (^) we have used a steepest descent estimate, 
in the limit R/Lq 0, in order to define the intermittent scaling exponents 
C{m). Intensity of intermittency depends on the departure of the C("^) ex- 
ponents from a linear behavior in m. 

In order to extend this description to the time domain, Procaccia and coll. 
have proposed to consider that two velocity fluctuations, both at scale R but 
separated by a time delay t, can be thought to be characterized by the same 
fragmentation process Wn^Loit) ~ W^i?,Lo(0) long as the time separation 
t is smaller then the "instantaneous" eddy-turn-over time of that scale, tr, 
while they must be almost uncorrelated for time larger then t^. Consider- 
ing that the eddy-turn-over time at scale R is itself a fluctuating quantity 
tr ~ R/{5vr) ~ R"^'^ we may write down |0: 

r / R\ ''(^^^^ / t \ 

C^'\R, R\t)^ j d^^,.,Sh) (-) F,,, [-J (5) 

where the time-dependency is hidden in the function f{x) which must be a 
smooth function of its argument (for example a decreasing exponential). 
From (^) is it straightforward to realize that at zero-time separation we 
recover the usual SF representation. It is much more interesting to notice 
that (^) is also in agreement with the constraints imposed by the non-linear 
part of the Navier-Stokes eqs. Indeed, to make short a long story (see 
for a rigorous discussion) we may say that under the only hypothesis that 
non-linear terms are dominated by local interactions in the Fourier space 
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we can safely assume that as far as power-law counting is concerned the 
inertial terms of Navier-Stokes equations for the velocity difference 6vr can 
be estimated to be of the form: 



dtSvRit) ~ O 



R 



(6) 



and therefore we may check that: 



dtCP''^{R,R\t)^ dfi,,,Sh) 



h{p+q) 



{^R 



P,Q 



t 

tr 



CP+^^'i{R,R\t) 
R 



(7) 



where of course in the last relation there is hidden the famous closure-problem 
of turbulence, now restated in term of the relation : ^Fp^q{t) ~ Fp+i Let 
us therefore stress that we are "not solving turbulence" but just building up 
a phenomenological framework where all the leading (and sub-leading, see 
below) scaling properties are consistent with the constraints imposed by the 
equation of motion^. 

In the following we shall show how the representation must be im- 
proved to encompass the most general multi-time multi-scale correlation 
{r,R\t). 



3 Single scale time correlations 

We shall first show in which respect the expression (|^) may not be con- 
sidered a satisfactory representation of single scale time correlation. The 
first comment that can be raised about (|^) is that it misses important 
sub-leading terms which may completely spoil the long-time scaling be- 
haviour: indeed, the main hypothesis that correlation C^''^{R, R\t) feels only 
the eddy-turn-over time of the scale R is too strong. It is actually cor- 
rect only when the correlation function has zero disconnected part, i.e. when 
lim^^oo ('^'^r(O)) {^'^Fti^)) = which is certainly false in the most general case. 
The problem is not only limited to the necessity of taking into account the 
asymptotic mismatch to zero given by the disconnected terms - which would 

^In order to really attack the NS equations in this framework one should dive into the 
structure of the .Fpq-functions in great detail: a problem which seems still to be far from 
convergence 1^. 
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be a trivial modification of (^) - because as soon as tlie disconnected part 
is present tlie whole hierarchy of fluctuating eddy-turn-over times from the 
shortest, t^, up to the largest, tig, must be felt by the correlation. 

Let us, for the sake of simplicity, introduce a hierarchical set of scales, 
In = 2"Lo with n = 0, . . . , rid, which span the whole inertial range such that 
h = ^rid = '^'"''^Lq, and let us simplify the notation by taking Lq = 1 and by 
writing m„ = 6vr in order to refer to a velocity fluctuation at scale r = l^- 
The picture which will allow us to generalize the time-multifractal represen- 
tation to the multi-time multi-scale case goes as follows. 
For time-delays, t ~ tm, typical of the eddy-turn-over time of the m-th scale 
we may safely say that the two velocity fluctuations follow the same frag- 
mentation process from the integral scale Lq down to scale m while they 
follow two uncorrelated processes from scale m down to the smallest scale 
in the game, n. In the multifractal language we must write that for time 
t = tm ± 0{tm) we have: 



«n(o) ~ i^;„(o)H^„,o(o)«o(o) ~ t (8) 

Unit) ~ W'l^{t)W^,Q{t) ~ P^;^(t)iy„,o(0)no(0) ~ f t (9) 

where the exponents h, h' , h" are independent outcomes of the same prob- 
ability distribution functions and where we have used the fact that in this 
time- window Wmflit) ~ const. 

Apart from subtle further-time dependencies (see below) we should therefore 
conclude that for time t ~ tm the correlation functions may be approximated 
as: 

C'nUtm) ~ (Wl^m) (W^S') (10) 

which must be considered the fusion-rules prediction for the time- dependent 
fragmentation process 0, ^, |^. Let us notice that this proposal has already 
been presented in and considered to express the leading term in the limit 
of large time delays tm oo; here we want to refine the proposal made in 
[|] showing that by adding the proper time-dependencies it is possible to 
obtain a coherent description of the correlation functions for all time-delays. 
The expression (|1^) summarizes the idea that for time delay larger than 
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but smaller then tm-i, velocity components with support on wavenumbers 
k < km-i did not have enough time to relax and therefore the local exponent, 
h, which describes fluctuations on those scales must be the same for both 
fields. On the other hands, components with support on wavenumbers k > 
km-i have already decorrelated for t > t^-i and therefore we must consider 
two independent scaling exponents h', h" for describing fluctuations on these 
scales. 

Adding up all this fluctuations, centered at different time-delays, we end 
with the following multifractal representation for C^'^(/„, ln\t) = C^'^()f:): 

n-1 „ / J \qhi / J \ph2 

m=l •' Xi'mJ \^mj 

Let US now spend a few words in order to motivate the previous expression. 
In the first row of ( [Til) we have explicitly separated the only contribution we 
would have in the case of vanishing disconnected part. This term remains 
the leading contribution in the static limit {t = 0) also when disconnected 
parts are non-zero. About the new terms controlling the behavior of the cor- 
relation functions for larger time we still have at this stage the most general 
dependency from all times entering in the game t, tn, tm- In practice one may 
guess a very simple functional form which is in agreement with all the above 
mentioned phenomenological requirements. In particular, we may simply as- 
sume that fp^q (^^, = fp^q (^-^^ with fp,q{x) being a function peaked at 
its argument x ~ 0(1) which must be exactly zero for x = and different 
from zero only in a interval of width 6x ~ 0(1). 

Let us now face the consistency of the representation (^) with the con- 
straint imposed by the equations of motion. By applying a time- derivative 
to a correlation C^'^(t) you produce a new correlation with by-definition 
zero-disconnected part, whose representation has thus no sub leading term 
(/p g = 0). When performing the time derivative on both sides of ([TT]) it 
is evident that - in order to accomplish the dynamical constraints - all the 
derivatives of the subleading terms must sum to a zero contribution. 
This is the first non-trivial result we have reached until now. If our represen- 
tation ( |lTl) is correct, we claim that all eddy-turn-over times must be present 
in the general single-scale correlation functions but strong cancellations of 
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all sub-leading terms must take place whenever disconnected contributions 
vanish. 

Let us finally notice that for time delays larger then the eddy-turn-over 
time of the integral scale we should add to the RHS of ([TT| ) the final expo- 
nential decay toward the full disconnected term (m^) (m^). 
In order to check this representation we have performed some numerical in- 
vestigation in a class of dynamical models of turbulence (shell models) 
Without entering in the details of this popular dynamical model for the tur- 
bulent energy cascade let us only say that within this modehng the approx- 
imation of local-interactions among velocity fiuctuations at different scales 
is exact and therefore no-sweeping effects are presents. This fact makes of 
shell models the ideal framework where non-trivial temporal properties can 
be investigated. 

In order to test the dependency of ( P^Tj ) from the whole set of eddy-turn- 
over times we show in Fig. 1 the correlation Cl^^{t) for two cases with and 
without disconnected part. As it is evident, the correlation with a non-zero 
disconnected part decays in a time-interval much longer then the characteris- 
tic time of the shell r„. This shows that it is not possible to associate a single 



time-scale Tn to the correlation functions of the form (11). In Fig. 2 we also 
compare the correlation when one of the two observable is a time-derivative 
with the correlation chosen such as to have the same dimensional properties 
but without being an exact time-derivative. Also in this case the difference 
is completely due to the absence (presence) of all sub-leading terms in the 
former (latter). 

3.1 Intermittent integral time-scales 

In the case when the disconnected part of the time correlation is absent, in 
the representation ( |TI1) all the subleading terms mutually cancel, leaving the 
fully-connected contribution alone. 

Under this condition and in presence of intermittency one expects anomalous 
scaling behavior for the integral time-scales, s^''^{R), characterizing the mean 
decorrelation time of fluctuations at scale R, deflned as [p!0| : 

^(..).^^ _ Iq- dtC^'%R,R\t) 
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exploiting the multifractal representation (|ri|) it is easy to show that: 



(R) - (j-) (13) 



where the exponents z(m) are fully determined in terms of the intermittent 
spatial scaling exponents: z{m) = 1 + ({m — 1) — C("^)- 
This prediction is in practice always very difficult to check: indeed full can- 
cellation of the subleading terms requires an extremely long time span, and 
since the cancellations affect dramatically the convergence of the time in- 
tegral, there is no chance of measuring with sufficient precision the z{m) 
exponents. 

In order to bypass this problem we devised an alternative way to extract the 
integral times. 

We introduce fluctuating decorrelation times at a scale R, deflned as the time 
interval t/"^'' in which the instantaneous value of the correlation has changed 
by a flxed factor A, i.e. in our octave notation: 

Un{ti)Un{ti + Ti) = X^^\Un{ti)\^. (14) 

At time t^+i = ti + Ti we repeat this procedure and we record the new 
decorrelation time Tj+i and so forth for an overall number of trials A^. The 
averaged decorrelation times can be thus deflned as 

= {T^lUrDe/iTlUnDe = {TlUnDt/ HUnDu 

where (• ■ ■)e stands for ensemble averaging over the trials and (■ ■ ■) rep- 
resents the usual time average]^. Since the multifractal description applies to 
time averages the averaged decorrelation times scale as r^™-* ~ with the 
same scaling exponents of the integral times s™. 

In Table 1 we report the observed numerical ({m) along with the observed 
and expected scaling exponents for t^"^\ showing a very good agreement. 



■^The relation between the e-average and the t-average is simply derived by observing 

that {KrT)t = iJ;[T\ur,rdt)/T:^iJ:,Tl\uniu)n/iJ:^T^ = {T^^njiT),. 
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4 Two-scales time correlations 



Let us now jump to the most general multi-scale multi-time correlation func- 
tions: 

C'^nit) = (KiO) ■ uUt)) (15) 

where from now on we will always suppose that describes the velocity 
fluctuation at the smallest of the two scales considered, i.e. N > n. It is 
clear that now we have to consider the joint statistics of two fields: the first, 
the slower, at large scale, Un{0) and the second, the faster, at small scale and 
at a time delay t, UNit)- 

Following the same reasonings as before we may safely assume that from 
zero time delays up to time delays of the order of the slower velocity field, 
t — tn, the velocity field at small scale feels the same transfer process of 
Un up to scale n and then from scale n to scale N an uncorrelated transfer 
mechanism: 

UN{t) = WN,n{t)Mt) ^iV,„(t)M„(0) - WN,n{t)KoiQH for < t < t„ 

(16) 

Similarly, for time delays within tn < tm < t < tm^i < to also the field at 
large scale n will start to see different transfer processes: 

MO) ~ W^:,to(0)«™(0) ~ K™(0)W^«,o(0)«o (17) 

UN(t) ~ WN,m{t)Um{t) ~ W^iV,mMm(0) ~ WN,rn{tW;^,o{0>0 (18) 

It is clear now, how we may write down the correlation for any time: 

m=l ^ ' 

where we want to stress that the sum in the above expression goes only up 
to the index of the largest scale n (see below) . In order to understand which 
would be a reasonable functional shape for the fp^q function we need to point 
out a few preliminary remarks. Once we have to cope with two-scale corre- 
lation functions is natural to suppose that the the time-delay, tnN = tn — tN, 
needed for a energy burst to travel from shell n to shell N will play an impor- 
tant role. Furthermore, as far as the time-decaying properties are concerned 
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we must require that only eddy-turn-over times from the slower time t„ up 
to the large-scale eddy-turn-over time to enter in the game. This is because 
only for time larger then tn the correlation is a true multi-time correlation. 
Indeed, for time-delay shorter then t„ only the field at small scale, un, is 
changing but always under the same large scale configuration, m„. 
The final /-shape may be therefore guessed as: fp^q (^-^, = fp^g ^ (*-*niv) 

Where again the function f{x) must be a function peaked for x ~ 1 and with 
a width Sx ~ 0(1). The matching of representation ([T9|) with the equa- 
tion of motion reveals some important dynamical properties. From simple 
time-differentiation we should have 



N 



(20) 



which seems to be in d2sagreement with the time-representation proposed 
([lyp because in the RHS of ( pUD does appear explicitly the fast eddy-turn- 
over time tN (through the dependency from l^)- Actually, the representation 
(|TU|) is still in agreement with the equation of motion because the dependency 
of (pOD from t^r is false: again, exact cancellations must take place in the RHS. 
The explanation goes as follows: in the multifractal language we may write 
UN(t) = W]\f^n(t)unit) and therefore 

j^Mt) = (^^^.-w) ""(^) + ^NAt) (21) 

but for time shorter then the eddy-turn-over, t„, of the large scale m„, the 
term W^At) (^) Un{t) is zero because the shell m„ did not move at all, while, 

once averaged, the first term of the RHS of ( pTj) becomes W^7v,n(^)) ) {un{t)) 
which also vanishes because of the total time derivative. The time deriva- 

five, dtC^'%{t) will therefore be a function which scales as instead of 

^•^^ — as simple power counting would predict. This may even be shown 
rigorously by evaluating the following averages: 

dt {<{t)uUt)) = = {{d,ul) O + « {d,u^^)) . (22) 

The previous exact relation forces one of the two correlation to not satisfy 
the simple multifractal ansatz because otherwise power-law counting would 
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be contradictory: 



mo o - + « {d^u\)) ^Sp^ (23) 

Now, in view of the previous discussion, we know that it is the correlation 
with the time derivative at small scale, {u'^{dtvF^)) , which does not satisfy 
the multifractal power law, but has the same scaling of {{dtuDu^j^) as our 
representation correctly reproduces. 

In order to test all these properties, we plot in Fig. 3 the typical multi- 
time multiscale velocity correlation C^'^„(t) for p = g = 1, = 6, = 6 — 13. 
As one can see the correlation has a peak which is in agreement with the delay 
predicted by (|1^), which saturates at the value TnN — Tn for N <^n. Let us 
also notice that due to the dynamical delay, tm,Ni the simultaneous multiscale 
correlation functions C^''„(0) do not show the fusion-rules prediction, i.e. 
pure power laws behaviors at all scales: 

// \^(^) 



Indeed, for t — > 0, the term m = n dominates in (|T9D (because it makes 
tnN minimum) and fp^q (— can be considered a constant only in the 
limit of large scale separation, n <^ N, while otherwise we will see finite-size 
corrections. 

The effect of the delay in multi-scale correlations is shown in Fig. 4 where we 
compare Cl}\{0) and Cl}\{TnN) (for N > n = 6) rescaled with the Fusion 
Rule prediction (^). The time delay TnN is the time of the maximum of 
CN^nit) computed from Fig. 3. We see that without delay, the prediction 
( ^4|) is recovered only for N ^ n with a scaling factor /i^i(— 1) ~ 0.83, while 
including the average delay TnN the Fusion Rule prediction is almost verified 
over all the inertial range. 

Of course the delay TnN is a fluctuating quantity and one should com- 
pute the average ([T9|) with fluctuating delays. In this case the dimensional 
estimate t^n = inU~^ — Inu]^ is somehow ill-defined, first of all being not 
positive definite. To find a correct definition for the fluctuating time delays 
is a subtle point which lays beyond the scope of the present Paper. 
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5 Conclusions 



In conclusion, we have proposed a multifractal-like representation for the 
multi-time multi-scale velocity correlation which should take into account all 
possible subtle time- dependencies and scale-dependencies. The proposal can 
be seen as a merging of the proposal made in - valid only for cases when 
the disconnected part is vanishing - and the proposal made in |p - valid only 
in the asymptotic regime of large time delays and large scale separation. 
Our proposal is phenomenologically realistic and consistent with the dynam- 
ical constraints imposed by the equation of motion. We have numerically 
tested our proposal within the framework of shell models for turbulence. 

A new way to measure intermittent integral-time scales, sP''^{R), has also 
been proposed and tested. 

Further tests on the true Navier-Stokes eqs. would be of first-order impor- 
tance. Furthermore, the building of a synthetic signals which would repro- 
duce the correct dynamical properties of the energy cascade would also be of 
primary importance [|ll|, . 

We thank R. Benzi, V. L'vov and I. Procaccia for useful discussions. 
This work has been partially supported by INFM (Progetto di Ricerca Avan- 
zata TURBO). G.B. thanks the "Istituto di Cosmogeofisica del CNR" , Torino, 
for hospitality. 
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TABLE CAPTIONS 

TABLE 1: Comparison between the integral time-scales intermittency 

exponents, z^, estimated from the measured spatial intermittent exponents, 
zl^^^ = 1 + ({m — 1) — C(?7T.), and from direct measuring via the "doubling- 
time" T, ^^"'"). 
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FIGURE CAPTIONS 



FIGURE 1: 

The time dependency of single-scale correlation functions, C„„(t), in two dif- 
ferent cases. The continuous line is the case with a non-zero disconnected 
part, Cnn{t) = (l'"n(0)| |Mn(t)|) — (|m„|)^, whilc the dashed line represents a 
case with vanishing disconnected part Cnn{t) = K((m„(0)m* (t))). Both corre- 
lations are rescaled to their value at t — 0. The scale is fixed in the middle of 
the inertial range, n — 12, and the eddy turn over time of the reference scale 
was Ti2 — 0.29. The average has been performed over approximately 500ri2, 
about 10 large eddy turn over times. The presence of subleading terms in 
Cnn is apparent. The remnant anticorrelation in Cnn, for t > T12 reveals a 
partial cancellation of subleading terms: full cancellation requires averaging 
over a time interval of many more large eddy turn over times. 

FIGURE 2: 

Comparison between the two correlations C„„(t) — kn{\un\'^{0)\un\^{t)), con- 
tinuous line. Dnn{t) — —{\un\'^{Q)^^^^{t)), dashed line. The two correla- 
tions have the same dimensional properties, but Dnn{t) decays faster due 
to cancellations of subleading terms. D„„ vanishes at zero delay because of 
stationarity and smoothness of the process Unit)- Scale and characteristic 
times as in Figure 1. 

FIGURE 3: 

Multi-time multi-scale correlation functions, Cn,N{t), for n = 6 and N = 
6, ■ ■ ■ , 13 (from bottom curve to top curve). Observe the saturation in the 
time-delays, Tn,N — Tn — tn ^ Tn — 0{1) when N increases. 



FIGURE 4: 

Lin-log plot of Multi-scale correlation C^f^i't) =< |wn(0)||wAr(i)| > rescaled 
with the Fusion Rule prediction: Cl'^N{t) /{Sj^S^/ S^) at fixed n = 6 and at 
changing N > n. The lower line represent the zero-delay correlation [t — 0), 
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the upper line is for the average delay t = Tq^n- 
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Ill 




,{iiim,) (.{III)) 

"'III y^'iii J 


1 


0.39 


-0.61 (-0.61) 


2 


0.72 


-0.68 (-0.67) 


3 


1.00 


-0.72 (-0.72) 


4 


1.26 


-0.75 (-0.74) 


5 


1.49 


-0.77 (-0.78) 


6 


1.71 


-0.78 (-0.78) 


7 


1.93 


-0.80 (-0.80) 


8 


2.13 


-0.80 (-0.80) 



Table 1: 
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